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Abstract. For the BCS equation with local two-body interaction W(x), we 
give a rigorous analysis of the asymptotic behavior of the critical temperature 
as A — > 0. We derive necessary and sufficient conditions on V(x) for the 
existence of a non-trivial solution for all values of A > 0. 



1. Introduction 

The BCS model has played a prominent role in condensed matter physics in the 
fifty years since its introduction [2] . Originally introduced as a model for electrons 
displaying superconductivity it has recently also been used to describe dilute cold 
gases of fermionic atoms in the case of weak interactions among the atoms [101 [HI 
HH [TJ [13l |H [3] . We will not be concerned here with a mathematical justification of 
the approximations leading to the BCS model, but rather with an investigation of 
its precise predictions. 

We consider the BCS equation for a Fermi gas at chemical potential and 
temperature T > 0, with local two-body interaction 2XV(x). Here, A > denotes 
a coupling constant, and the factor 2 is introduced for convenience. Because of 
the many different applications of the BCS equation, it is important to keep the 
discussion as general as possible. Our only assumption on the interaction potential 
V will be that it is real-valued and V € £ X (R 3 ) H L 3 / 2 (R 3 ). 

It was shown in [5] that the existence of a non-trivial solution to the BCS gap 
equation 

with E(p) — y/ (p 2 — /i) 2 + |A(p)| 2 at some temperature T > is equivalent to 
the fact that a certain linear operator has a negative eigenvalue. Here, V(p) = 
(27r)~ 3 / 2 J V(x)e~ lpx dx denotes the Fourier transform of V(x). In particular, it 
was shown that this property holds for T less than a certain critical temperature, 
which we denote by T c , whereas there are no non-trivial solutions to Eq. (fTJ) for 
T > T c . According to the usual interpretation of solutions to the BCS gap equation, 
the system displays superfluid behavior for all temperatures T < T c , while it is in 
a normal phase for T >T C . The analysis in [5] shows that T c is non-zero for purely 
attractive (i.e., non-positive) V, and exponentially small in A. 
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The fact that the critical temperature in the non-linear BCS equation can be 
expressed in terms of spectral properties of a linear operator allows for a more 
thorough investigation of its properties. This is the purpose of this paper. In 
particular, we shall be concerned here with the asymptotic behavior of T c at weak 
coupling, i.e., for small A. We shall derive necessary and sufficient conditions on V 
for the positivity of T c for all A > 0, as well as its precise asymptotics as A — > 0. 
The precise statement of our results is given in Theorem Q] below. 

The linear operator one is led to analyze is of the form Kt,^ + XV where Kt,^ is 
a multiplication operator in momentum space that represents an 'effective' kinetic 
energy. By a modification of the Birman-Schwinger principle we need to study the 
diverging part of the compact operator 



as T — > 0. Note that, if V is not of definite sign, then the latter operator is not self- 
adjoint and standard perturbation arguments based on the variational principle will 
fail. Still we are able to give a variational characterization for the leading behavior 
of the critical temperature in the weak coupling limit. 

Our analysis is somewhat similar in spirit to that of the lowest eigenvalue of 
the Schrodinger operator p 2 + XV in two space dimensions, see |15j . This latter 
case is considerably simpler, however, since p 2 has a unique minimum at p = 0, 
whereas Kt^{p) takes its minimal value on the Fermi sphere p 2 = /i. Technically, 
this is reflected in the fact that the singular part of the Birman-Schwinger operator 
(sgn V)|y| 1/2 (p 2 + 2") _1 |V| 1/2 is of rank one in contrast to that of ©, which is of 
infinite rank. In particular, the difficulties stemming from the non-selfadjointness 
are not present in the case of p 2 + XV. 

We would like to emphasize that our approach is not restricted to the kinetic 
energy Kt,^ appearing in the BCS model, but can be adopted to any symbol 
vanishing on a manifold of codimension one or higher. Operators of this form arise 
naturally in various fields of Mathematical Physics, e.g. in the quantum-mechanical 
description of particles in a homogeneous magnetic field or in the analysis of trapped 
modes in elasticity theory [5j [16] . 



According to the analysis in [8], the critical temperature in the BCS model is, 
in appropriate units, given by the following expression. 

DEFINITION 1. For /x € E and T > 0, let K T ^ be the multiplication operator 
in momentum space 



Let V G L^R 3 ) n L^M 3 ) be real-valued. The critical temperature in the BCS 
model is given by 



(sgn^)|y|V^-i|y|V2 



(2) 



2. Mail results and discussion 




T C (V) = inf {T > : inf spec {K T ,„ + V) > 0} . 



(3) 



More precisely, it was shown in [5] that Eq. (fTJ) has a non-trivial solution for 
T < T C (XV), whereas for T > T C (XV) it doesn't. Note that K T ^ > 2T, and that 
the essential spectrum of Kt^ + V is [2T, oo). Hence, in case T C (V) > 0, it is the 
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largest T such that Kt,^ + V has a zero eigenvalue. Note also that Kt,h becomes 
\p 2 - fj,\ as T -> 0. 

We assume that /i > henceforth. For weak potentials V, the critical temper- 
ature is determined by the behavior of the potential on the Fermi sphere f2 M , the 
sphere in momentum space with radius y/JL. We denote the Lebesgue measure on 
fi^ by dw. 

Let Vfi : L 2 (f2 M ) — > L 2 (f2 M ) be the self-adjoint operator 

CM (P) = tA^A= I V(p - q)u(q) dw{l) ■ (4) 



We note that V M is non-vanishing if V(p) does not vanish identically for \p\ < 
Since V G L 1 (R 3 ) by assumption, V is a bounded continuous function, and hence 
V M is a Hilbert-Schmidt operator. It is, in fact, trace class, as will be shown below, 
and its trace equals ^§ / V(x)dx. 

Let a^{V) — inf spec (V M ) denote the infimum of the spectrum of V M . Since V M is 
compact, we have a M (V) < 0. Note that, in particular, a^{V) is negative if the trace 
of V M is negative, that is, a M {V) < if ^(0) = {2tt)- 3 ^ J V(x)dx < 0. Moreover, 
by considering a trial function that is supported on two small sets on the Fermi 
sphere separated a distance \p\, it is easy to see that a^V) < if |V"(p)| > V"(0) 
for some p with \p\ < 2y/Ji. 

Our main result concerning the critical temperature ([3]) is as follows. 

THEOREM 1. Let V G £ 3/2 (R 3 ) n L^R 3 ) be real-valued, and let A > 0. 

(i) Assume that a M (V) < 0. Then T C (XV) is non-zero for all A > 0, and 

lim A In — /f — - = ; — - . (5) 

a^o T C (\V) a„{V) V ; 

(ii) Assume that a^(V) = 0. IfT c (XV) is non-zero, then hi(n/T c (XV)) > cA~ 2 
for some c > and small A. 

(hi) If there exists an e > suc/i that a^(V — e\V\) — 0, then T C (\V) — for 
small enough A. 

Note that Eq. ([5]) implies that, in case a^{V) < 0, the critical temperature has 
the asymptotic behavior 

T C (XV) ~ ^e 1 l( Xa ^ v )) 

in the limit of small A. On the other hand, if a^(V) = then part (ii) of Theorem[T] 
implies that T C (XV) is at most as big as e - const / A for some positive constant. If 
a^V) remains zero if e|V(x)| is subtracted from V(x), then T C (XV) = for small 
enough A, and there is no superfluid phase at weak coupling. 

Although we restrict our attention to local potentials V here, we remark that a 
similar analysis can be applied in the case of non-local potentials as well. 

2.1. Radial Potentials. In the special case of radial potentials V, depending only 
on | a; |, the spectrum of V M can be determined more explicitly. Since com- 
mutes with rotations in this case, all its eigenfunctions are given by spherical har- 
monics. For £ a non-negative integer, the eigenvalues of V M are then given by 
J V(x)\je(yfjl\x\)\ 2 dx, with jg denoting the spherical Bessel functions. These 
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eigenvalues are (2£+ 1) fold degenerate. In particular, we then have 

a»(V) = inf ^§ J V(x) \MVt\x\)\ 2 dx ■ 

in the case of radial potentials V. We remark that J2 ee f S (2£ + l)li«( r )| 2 = 1> hence 
the expression for the trace of V M stated above is recovered. 

If V is non-positive, it is easy to see that the infimum is attained at I = 0. 
This follows since the lowest eigenfunction can be chosen non-negative in this case, 
and is thus not orthogonal to the constant function. Since jo(r) — sin(r)/r, this 

means that a^V) — (2ir 2 y/JI)^ 1 J V(x) sm WS^- dx for radial potentials V with 
non-positive Fourier transform. 

In the limit of small we can use the asymptotics je(r) w r e /(2£ + 1)!! to 
observe that, in case fV(x)dx < 0, a^V) ~ fV(x)dx as /i — > 0. Note that 
(X/Att) JV(x)dx is the first Born approximation to the scattering length of 2XV, 
which we denote by ag. Thus, replacing Xa^(V) by 2y / /7ao/ - 7r and writing // = kf, 
we arrive at the expression T c ~ e n '( 2kfa °> for the critical temperature, which is 
well established in the physics literature [12l [3] . 

In the remainder of this paper, we shall give the proof of Theorem [1] 

3. Proof of Theorem [1] 

Note that in case T c > 0, the essential spectrum of Kx c ,p, + XV starts at 2T C > 0, 
and hence T c is the largest T such that is an eigenvalue of Kt,^ + XV in this case. 
Therefore there exists an eigenstate £ L 2 (R 3 ) such that K^.^ip} = —XV\ip). 
For a (not necessarily sign-definite) potential V(x) let us use the notation 

V(x) 1/2 = {sgnV(x))\V(x)\ 1/2 . 

The Birman-Schwinger principle then implies that \<p) = V 1 ' 2 ^) satisfies Bt c \<p) — 
— \ip), where 

B T = XV 1/2 K~^\V\ 1/2 . (6) 

Conversely, if B T \ip) = -\tp) and = K T ^\V\^ 2 \ip) , then \ip) e L 2 {R 3 ) and 
^T,ii\ip) — —XV\ip). The existence of a zero eigenvalue for K T ^ + XV is thus 
equivalent to the fact that Bt has an eigenvalue —1. Note that Bt is not a self- 
adjoint operator, however. 

With the aid of the Birman-Schwinger operator Bt, we can thus state the fol- 
lowing alternative characterization of the critical temperature T C (XV). 

Lemma 1. For any T > 0, the Birman-Schwinger operator Bt defined in (0|) is 
Hilbert- Schmidt and has real spectrum. If T C {XV) > 0, the smallest eigenvalue of 
Bt c equals —1. Moreover, in case T C (XV) = 0, the spectrum of Bt is contained in 
(—1, oo ) for any T > 0. 

Proof. The Hilbert-Schmidt property follows from the Hardy-Littlewood-Sobolev 
inequality [TO Thm. 4.3], using that V S L 3/2 (R 3 ) and that K T ^ > const (1 +p 2 ). 
Moreover, Bt is the product of a self-adjoint operator (multiplication by sgn(V(a;))) 
and a non-negative operator, hence it has real spectrum. 

We have already shown above that —1 is an eigenvalue of Bt c in case T c > 0. 
Moreover, because of strict monotonicity of Kt.^ in T, —1 is not an eigenvalue of Bt 
for all T > T c . This implies that Bt c has no eigenvalue less than —1, for otherwise 
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there would be a T > T c for which Bt has eigenvalue —1 since the eigenvalues of 
Bt depend continuously on T and approach as T — > oo. 

In the same way, one argues that Bt does not have an eigenvalue less than or 
equal to -1 if T c = 0. □ 

Let J be the unitary operator that multiplies by sgn(V(x)). To be precise, we 
define sgn(V(x)) = 1 in case V(x) = 0. Moreover, let X denote the self-adjoint 
operator on L 2 (R 3 ) with integral kernel 

x^y) = \v^ ^f x - v \ v{y)\^- 

\x — y\ 

We note that the X is a non- negative trace-class operator, with trace tr[X] = 
70s J\V(x)\dx. Hence also JX is trace-class, and trfJX] = ^ JV(x)dx. Define 
F T by 

B T = Aln(l + ^) JX + XY T . (7) 

We have 

Lemma 2. Let V G L 1 {W 3 ) n L 3 / 2 (R 3 ). Then, for any T > 0, the operator Y T 
defined in ^ is Rilbert- Schmidt, and its Hilbert- Schmidt norm is bounded uniformly 
in T, i.e., sup T>0 tr[l^tlr] < oo. 

The proof of this lemma will be given in the next section. Lemma [2] shows that 
the singular part of the operator Bt as T — * is given by JX. This observation 
will enable us to recover the exact asymptotics of T C (XV) as A — * 0. 

The operator JX is closely related to V M defined in In fact, the two operators 
are isospectral. 

Lemma 3. The spectrum of JX on L 2 (R 3 ) equals the spectrum ofV^ on L 2 (f^). 

Proof. Let A : L 2 (R 3 ) i-> L 2 (VL^) denote the operator which maps ip G L 2 (K 3 ) to 
the Fourier transform of |y| 1 / 2 -0, restricted to the sphere Note that | 1 1//2 -0 G 
L 1 (M. 3 ) and hence it has a bounded and continuous Fourier transform. Moreover, 
let B : L 2 (f^) h-> L 2 (M 3 ) be defined by 

(Bu)(x) = V(x)^j^-LJ^ u(pyv* MP) ■ 

Using the fact that J n e tpx duj(p) — 47r v //7|a;| _1 sin ^/yu|a;| it is easy to see that 
JX = BA, while AB = V M . Hence they have the same spectrum, except possibly 
at zero. Indeed, if AB\f) = X\f) with A ^ 0, then \g) = B\f) ^ and BA\g) = X\g). 
Since both operators are Hilbert-Schmidt operators on infinite-dimensional spaces, 
is an element of both spectra. □ 

We now study the behavior of the spectrum of JX under small perturbations. 
We will show that for a > the spectrum of 

aJX + XY T (8) 



differs from the spectrum of aJX by at most 6(V aA) + 6(A), uniformly in T. 
Here and in the following, we use the notation 0(i) to indicate an expression that 
is bounded as ct < 0(t) < Ct for constants < c < C. 
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Pick a z that stays away a distance d from the spectrum of aJX . By expanding 
in a Neumann series, we see that aJX + XYt — z has a bounded inverse provided 



X\\Y T \ 

We have 



1 



aJX 



< 1. 



aJX-z z z aX 1 / 2 JX 1 / 



2 



Z 



Since X l / 2 JX l l 2 is a self-adjoint operator having the same spectrum as JX, we 
can bound (aX 1 / 2 J X 1 / 2 — z)\\ < 1/d for any z a distance d away from the 
spectrum of a JX. We conclude that || (a JX — 1| < 1/d + a||X|j/d 2 . Hence z 
is not in the spectrum of © if d > 0(-\/Aa) + 9(A). 

Since the spectrum of aJX + XYt depends continuously on A, we have thus 
proved the claim. In particular, it follows that the lowest eigenvalue of {8} equals 
the lowest eigenvalue of aJX plus terms that are at most of order O(vaA) + 0(A). 

We now have the necessary prerequisites to give the proof of Theorem [TJ 

Proof of Part (i). According to LemmaEl we have a^(V) = inf spec JX. Assume 
now that inf spec JX < 0. Since Yt is bounded uniformly in T, we see that the 
spectrum of Bt = Aln(l + /x/2T) JX+ XYt becomes arbitrarily negative for T — ► 0, 
and hence T C (XV) > for any A > 0. Moreover, Aln(l + /i/2T c ) is bounded away 
from zero as A — > 0. 

We have shown above that the lowest eigenvalue of Bt is bounded from above 
and below by Xa^V) ln(l + /x/2T) + 9(\/A), uniformly in T. Since at T = T c this 
lowest eigenvalues equals — 1, we conclude Eq. ([5|). □ 

Proof of Part (ii). For T > 0, let a — Aln(l + n/2T) for simplicity. Under the 
assumption that the spectrum of JX is non-negative, the lowest eigenvalue of Bt — 
aJX + XYt is bigger than — 9(v / aA), as shown above. This immediately implies 
that Bt can only have an eigenvalue —1 if aX > 6(1), or ln(/i/T) > 0(1/A 2 ) for 
small A. □ 

Proof of Part (Hi). Let again a = Aln(l+/z/2T), and recall that Bt — aJX + XYr- 
Since the operator 1 + XYt is invertible for small enough A, we are able to rewrite 

1 + B T = (1 + XY T ) (l + (1 + XY T )~ 1 aJX) . 

Hence 1 + Bt does not have a zero eigenvalue for any a > if the spectrum of 
(1 + XYt)^ 1 J X is non-negative. Note that JYt is self-adjoint, since JBt is self- 
adjoint and J 2 = 1. Hence (l + XYr) -1 JX has the same spectrum as the self-adjoint 
operator 

This operator is non-negative for small A if X 1 / 2 JX 1 / 2 > eX for some e > 0, since 
then 

© = X 1 I 2 JX 1 ' 2 - XX 1/2 Y T = \ 1V X 1 ' 2 >X(e- X\\Y T \\ ||(1 + AFt) _1 ||) . 

J + A j Yt 

Note that the range of X is dense in the range of X 1 / 2 , and hence it is enough 
to check the inequality J > e on the range of X. Let be in the range of X, 
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i.e., = X\4>) for some \<f>) e L 2 (R 3 ). Then (ip\Jip) > e(ip\ip) is equivalent to the 
statement that, for |x) = |U| 1/2 |0), 

/ x(p)V(p - q)x{q) duj(p) dw{q) 

J 11 n X 11 fj, 

>e/ x(p)\V\(p - q)x(q) dui(p) du{q) . 

This, in turn, is equivalent to a^iV — e\V\) = 0. Under this assumption, we have 
thus shown that, for small enough A, the operator Bt does not have an eigenvalue 
— 1, for arbitrary T > 0. Together with Lemma [TJ this proves the claim. □ 



4. Proof of Lemma 

By scaling we may assume that fj, = 1, and we set Kt,i = Kt for simplicity. 
The operator Kt can be rewritten as 

K T (p) = (\p 2 -l\+2T)g(\p 2 -l\/T) 
where git) = t(l + e _i )/((t + 2)(1 — e _t )). The integral kernel of K^ 1 is given by 

I r e ip[x-y) 
K T^V) =J2^J M3 (\p2-l\ + 2T)g(\p2-l\/T) dp - 

We decompose K T {p) as K^p)- 1 = L^(p) + M T 1] {p), where L^\p) = 9{V2 - 
\p\)K T (p)~ l and M { T 1] [p) = 9{\p\ - \f2)K T {p)~ x . Since b = w£ t g(t) > one has 

M T l \p) < 9{\p\ - sfyb-^p 2 - . 

Using that V € L 3 / 2 (R 3 ), we find with the aid of the Hardy-Littlewood-Sobolev 
inequality pH Thm. 4.3] that ||U 1 / 2 a4 1) |U| 1/2 ||2 is bounded independently of T. 
Here, || • || 2 = (tr| • | 2 ) 1/2 denotes the Hilbert-Schmidt norm. 
Note that the integral kernel of iSp is given by 

1 f^ 2 k smk\x — y\ 

1 ■ ak . 



2^J (\h?-l\+2T)g(\k*-l\/T) \x-y\ 
We further decompose Lj) = L^) + M T 2 \ where 

r (2)/ ^_ 1 f k wak\x-y\ 

L T^y)- 2n 2j o { \ k 2_ ll+2T) \ x _ y \ dk - 

Estimating | sin/c|a; — y\\ < v2|x — y\ for k < \[2 and changing variables one easily 
finds 

This is bounded independently of T since l/g(i) — 1 ~ 2/t as t — > oo. Since 
V S L 1 (R 3 ), we can bound HV^M^lVl^lla < f\V(x)\dx sup x>1/ \M%\x, y)\, 
and hence we see that also \\V 1/2 M T 2) \V\ 1/2 \\ 2 is bounded uniformly in T. 
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(2) (3) (3) 

Finally, we decompose L T = L T + M T , where 

{X ' y) -2«iJ (|P-l| + 2T) dfc \x-y\ 

Since | sin a — sin6| < \a — b\ one easily sees that 

rV2 



Again, since V G L^IR 3 ), || V 1/2 M^ ] \V\ 1/2 \\ 2 is uniformly bounded. This completes 
the proof. □ 
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